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Parton Distributions for LHC 
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I review the current status of the nucleon PDFs determined from global fits with emphasis on the impact of 
recent experimental data and the remaining theoretical challenges. 
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The detailed knowledge of the momentum dis- 
tributions of partons in the nucleon is quite nec- 
essary for phenomenology of the hard scatter- 
ing processes. Since the time, when the par- 
ton model was established the parton distribu- 
tion functions (PDFs) were permanently in the 
focus of the QCD studies. With a tremendous up- 
grade of the Tevatron collider luminosity and the 
start-up of the LHC at a record collision energy 
we need, however, an even better understand- 
ing of the PDF details in order to meet the im- 
proving accuracy of the experimental data. The 
PDFs are extracted from the hard-scattering data 
with the QCD evolution taken into account. At 
the present accuracy of the experimental data, 
moreover for the foreseen precision at the LHC, 
the QCD corrections at least up to the next-to- 
leading order (NLO) are required. Furthermore, 
the next-to-next-to-leading order (NNLO) correc- 
tions are necessary for many important processes, 
as W, Z-bosons, Higgs boson and top-quark pro- 
duction. 

The choice of the data used to constrain the 
PDFs is often defined by the accuracy of the avail- 
able theoretical calculations. For example, in or- 
der to include the DIS data with low momentum 
transfers Q 2 into the PDF fit, one has to take 
into account the higher-order QCD corrections 
and to take care about the higher-twist terms for 
the small-mass final state kinematics. Likewise, 
modifications of the nucleon wave function in nu- 
clei have to be taken into account once nuclear 
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data are used in the analysis. 

There are several nucleon PDF sets, which were 
kept updated during the recent years, cf. Table[TJ 
All these PDF sets are based on the inclusive 
deep-inelastic-scattering (DIS) data collected in 
fixed-target experiments and the data from the 
HERA collider. The HERA data mostly provide 
constraints on the gluon and sea quarks distri- 
butions at small values of the Bjorken variable x, 
which are particularly important for the interpre- 
tation of the future data of the first LHC run. Re- 
cently the inclusive data obtained by the HI and 
ZEUS experiments of Run I of the HERA collider 
were combined into the common data set j6] . The 
advantage of such an approach is mutual cross- 
calibration of both experiments, which leads to 
an essential improvement in the data accuracy. 
The slope of the updated HERA data on the in- 
clusive structure function (SF) Fi with respect to 
ln(Q 2 ), obtained from a model-independent fit, is 
given in Fig. [TJ It is in good agreement with 
the NNLO predictions based on the ABKM09 
PDFs [7 . However, the constant term obtained 
in the same fit systematically overshoots the pre- 
dictions. The combined HERA data in general lie 
above the separate data sets of Ref. [8] used in the 
PDF fits before. This change in the data requires 
a modification of the small- a; shape of the PDFs. 
In particular, for the version of the ABKM09 fit 
based on the combined HERA data, the gluon 
distribution at small x is shifted to lower values, 
cf. Fig. [2j In the updated version of the NLO 
fit by the CTEQ collaboration a decrease of the 
gluon distribution at small x is also observed [2]. 
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Table 1 

Recently published nucleon PDFs sets with a brief description of the data used in the fit, the theoretical 
accuracy and the factorization scheme employed to model the heavy-quark DIS contribution. 



Name Data used QCD approximation Scheme Reference 



MSTW 


DIS+DY+jets 


NLO/NNLO 


GMVFN 


DP 


CTEQ 


DIS+DY+jets 


NLO 


GMVFN 




NNPDF 


DIS+DY+jets 


NLO 


ZMVFN 


3 


JR 


DIS+DY 


NLO/NNLO 


FFN 


A 


ABKM 


DIS+DY 


NLO/NNLO 


FFN 




HERAPDF 


DIS 


NLO 


GMVFN 


m 



However, the MSTW group found much smaller 
variations of the PDFs jlOj . 

The PDF sets of Refs. [TT2"|5] were fitted also 
to the data on the and Z production from 
hadron colliders. These collider data provide im- 
portant constraints on the PDFs at large val- 
ues of the factorization scale. The W^, Z rates 
and distributions are calculated at NNLO |ll|12j 
and are considered as one of the LHC luminosity 
monitors. However, the most recent and accu- 
rate Run II Tevatron data 13 on the W charge 
asymmetry and the related data on the lepton 
charge asymmetry are not in a good agreement 
with the NNLO calculations of Ref. [12] based on 
the different PDFs, particularly for the case of 
the electron charge asymmetry data. The best 
agreement was observed for the case of the JR 
PDFs of Ref. 0], cf. Fig. [3] The MSTW and 
ABKM09 predictions overshoot the data in gen- 
eral. The NNPDF collaboration also obtained 
a poor description of the I-F-asymmetry data in 
their NLO fit 3 . Due to the QCD evolution of 
the PDFs the kinematics of the W-production 
at the Tevatron is constrained by the fixed tar- 
get DIS data at x = 0.1 -r- 0.2. Therefore the 
PF-asymmetry predictions are mostly sensitive to 
the valence quark isospin asymmetry, (u — d). 
This distribution can be constrained by the data 
on isospin-asymmetric combination of the proton 
and the neutron DIS SFs. However, the neutron 
target is available in a form of a bound state only. 
Therefore the neutron SFs unfolding requires to 
take into account nuclear effect. The DIS nuclear 
corrections are non- negligible even for the lightest 



nuclei, as the deuteron. Furthermore, they can- 
not be fully calculated using present theoretical 
methods in nuclear physics and are often fitted 
to nuclear data instead, cf. e.g. Ref [13]. The 
MSTW group included into their fit the deuteron 
correction in a form of a mo del- independent func- 
tion and fitted it to the global data set pi)]. In 
this way they found marginal improvement for 
the charged-lepton asymmetry. Moreover, a rea- 
sonable nuclear model cannot justify the form of 
the deuteron correction obtained. Besides, the 
data on the lepton asymmetry obtained by two 
Fermilab experiments, CDF and DO, are not in 
agreement. In view of these problems the CTEQ 
group suggests two variants of their recent PDF 
release, with and without Run II Tevatron data 
included [2]- 

The Tevatron jet inclusive production data 
were also used in Refs. |l|2|3j . Historically, these 
data were always considered as a unique source 
of the information about the gluon distribution 
at large x. Since the NNLO corrections to the 
hadronic jet production are not yet known these 
data can be used in a fully consistent way for ex- 
traction of the NLO PDFs only. With the Run I 
data of Ref. [15] included into the PDF fit, the 
CTEQ and MSTW groups observed a substan- 
tial increase of the large- a; gluons [16 17 . In fact, 
this effect was a manifestation of the tension be- 
tween the Run I jet data and the DIS data at 
large x. The accuracy of the data obtained in the 
Run II experiments has been greatly improved, in 
particular due to the better jet energy calibration. 
Furthermore, these data lie lower than those from 
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Figure 1. The slope di^/dlnQ 2 (upper panel) 
and the constant term (lower panel) for the com- 
bined HERA Run I data on the inclusive SF F2 
in comparison to the predictions of the ABKM09 
fit [6] . The impact of the correction on the contri- 
bution of the longitudinal SF F^, which were em- 
ployed to extract the values of F2 from the cross 
sections data is given by stars. Figure taken from 
Rcf. [Jj. 



Run I and are in much better agreement with the 
other data used in the PDF fits. As a result, the 
impact of the Run II results on the PDFs is re- 
duced, cf. Fig. 01 Due to this change of trend, 
the value of the strong coupling constant a s ex- 
tracted from the Tevatron data is shifted down- 
ward. The shift in the value of a s is essential 
since in leading order QCD the jet production 
cross section is 0(a 2 ). The value of a s (Mz) = 



0.1161 



+0.0011 

-0.0018 



was obtained in the recent analysis 
of the DO data performed at the NLO with ac- 
count of NNLO threshold corrections [18J . Tak- 
ing into account that the a^ NLO is smaller than 



a 



NLO 



this result is in a good agreement with 
the value of af NLO (M z ) = 0.1147(12) obtained 
in the updated version of the ABKM fit [7J and 



withaf ALU {M Z ) = 0.1141 



+0.0020 



obtained in the 



Figure 2. The la bands for the NNLO gluon dis- 
tribution obtained in the updated version of our 
fit (blue) compared to the same for ABKM09 fit 
(black) and the version of the updated fit with- 
out partial NNLO corrections to the heavy-quark 
electro-production of Ref. [9] taken into account 
(red) at the factorization scale of 9 GeV 2 . Figure 
taken from Ref. [7j. 



analysis of the non-singlet DIS data [T5]. The 
value of of LO (M z ) = 0.1215 obtained in the up- 
dated version of the MSTW fit [TU] is much larger 
than both these determinations. For comparison, 
the value of a s extracted from the e + e~ data is 
~ 0.119 [20]. However, it is worth to note that in 
the elaborated analysis of the global e + e~ data on 
the thrust distribution performed in the NNLO 
approximation with account of the power correc- 
tions to fragmentation the value of a s (Mz) = 
0.1135± (0.0002) exp . ± (0.0005) hadr . ± (0.0009) pcrt . 
was obtained [21] . This is much lower than the 
average of Ref. [20] and in good agreement with 
the determinations based on the DIS data. 

The HI and ZEUS data on the DIS charm pro- 
duction obtained in Run I of the HERA collider 
have also been combined into a common data 
set [33] recently. Heavy-quark electro-production 
is particularly sensitive to the shape of the gluon 
distribution at small x. Therefore it gives a 
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Figure 3. The Run II electron charge asymme- 
try obtained by the DO collaboration in compari- 
son with the NNLO predictions based on the JR 
PDFs [4]. Figure taken from Ref. [12]. 



complimentary constraint for the global PDF 
fits. Moreover, the detailed understanding of 
the heavy-quark DIS electro-production mecha- 
nism is necessary for the small- x inclusive SF de- 
scription. In the global PDF fits the zero-mass 
variable-flavor-number (ZMVFN) scheme is often 
used to model the heavy-quark contribution. In 
this scheme all quarks are assumed to be massless, 
which greatly simplifies the calculation. How- 
ever, the low-Q 2 part of the DIS SFs cannot be 
calculated within this concept due to the power 
corrections and different constant terms occur- 
ring in the coefficient functions. The ZMVFN 
scheme commonly overestimates the heavy-quark 
contribution at small Q 2 . Therefore it cannot be 
employed for the analysis of realistic DIS data. 
Instead of this the general-mass variable-flavor- 
number (GMVFN) scheme is used. The GMVFN 
schemes include a peculiar modeling of the low- 
Q 2 DIS region. This allows to overcome prob- 
lems of the ZMVFN scheme. However, since the 
modeling cannot be performed on solid theoret- 
ical grounds the GMVFN scheme cannot be de- 
fined in a unique way. Various variants of the 



Figure 4. The large- a; gluon distribution obtained 
in different variants of the NLO MSTW08 fit. 
Figure taken from Ref. pQ. 



GMVFN scheme were considered, cf. Ref. [23 a 
and references therein. Further modifications of 
the ACOT scheme [23] and the updated version 
of Thome's prescription [25] were suggested. The 
latter assumes a very flexible form of the semi- 
inclusive heavy-quark SF F£ at low Q 2 , which is 
provided by several additional parameters. An 
additional flexibility of model [35] would evi- 
dently allow better description of the data within 
this approach. However, additional parameters 
introduced in the GMVFN schemes to model the 
DIS SFs at low Q 2 have to keep the factoriza- 
tion scheme consistency. Otherwise the PDFs ob- 
tained in the fit based on the GMVFN scheme 
might correspond to some peculiar factorization 
scheme, but not to the MS scheme, commonly 
used in the QCD calculations, cf. Ref. [26 . In 
contrast, the BMSN prescription of the GMVFN 
scheme [37J gives a consistent MS description of 
the heavy-quark contribution in the transition re- 
gion. At the same time it provides a smooth tran- 
sition between the FFN scheme at low Q 2 and the 
ZMVFN scheme at large Q 2 , cf. Fig. [5J These 
advantages of the BMSN ansatz were also em- 
ployed in the FONLL approach to DIS considered 
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Figure 5. Comparison of F£ in different schemes 
to Hf- and ZEUS-data. Solid line: GMVFN 
scheme in the BMSN prescription, dash-dotted 
line: FFN scheme, dash line: ZMVFN scheme. 
Figure taken from Ref. [5j. 



in Ref. [28]. 

In the fixed-flavor-number (FFN) scheme the 
heavy-quark mass dependence follows from a 
complete fixed-order computation without any 
modeling even at low Q 2 . The coefficient func- 
tions of the heavy-quark semi-inclusive SFs are 
calculated in complete form up to the 0(a 2 ) 
(NLO) semianalytically. At the values of 
Q 2 3> the fixed-order FFN expressions are 
insufficient since they do not include the big- 
log terms of 0(a™ ln m (<5 2 )), which might be 
important even for large n. In the ZMVFN 
scheme these terms are resummed in a natural 
way through the QCD evolution of the heavy- 
quark PDFs. This is often considered as an ad- 
vantage of the ZMVFN schemes as compared to 
the FFN one. Nonetheless, these terms are par- 
tially taken into account in the NLO expressions 
of Ref. [35] and the numerical effect of the re- 
maining big-log terms is marginal for the realistic 
data kinematics. This observation was first made 
in Ref. [30] and later confirmed in Ref. [5] for the 



Figure 6. Comparison of PDFs obtained in 
Ref. [5J for the GMVFN scheme of Ref. [25] and in 
the variant of this fit for the FFN scheme. Figure 
taken from Ref. 1221. 



case of recent data on F£ and Ff accumulated 
by the HERA experiments. In particular, the dif- 
ference between the PDFs obtained in the FFN 
fit and in the GMVFN fit based on the BMSN 
prescription was found to be marginal [SJ. For 
the PDF fit, which includes the combined HERA 
data on F£ , the FFN scheme and the different 
variants of the GMVFN schemes provide equally 
good description of the data, cf. Ref. [22] . At the 
same time the gluon distributions obtained for the 
GMVFN scheme of Ref. [25] and the FFN PDFs 
are quite different, particularly for the small- x 
gluons, cf. FigM The ZMVFN gluons at small x 
must take lower values than for the FFN scheme 
by the scheme definition. However, this does not 
fully explain the effect given in Fig. [6] Meanwhile 
this effect might be related to the discrepancy 
in the gluon distributions obtained in the fits of 
Ref. [I] and Ref. [5j. 

The gluon-gluon collision is the dominant 
mechanism for the production of Higgs bosons 
with masses 100 < Mh < 300 GeV in hadron- 
hadron collisions. Furthermore, the leading or- 
der gluon-gluon term is 0(a 2 ). This makes the 
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Figure 7. The la error band for the Higgs-boson 
production cross sections at Tevatron and the 
LHC employing the ABKM PDFs [5] in compari- 
son with the central values for the case of MSTW 
PDFs |T] (dash-dotted lines). Figure taken from 
Ref. 0. 



Higgs production rate very sensitive to the vari- 
ation of the gluon shape and the value of a s . 
The Higgs production rates at the Tevatron col- 
lider and at LHC energies calculated with account 
of the NNLO corrections of Ref. [3T] are given 
in Fig. [7J The estimates based on the NNLO 
ABKM09 PDFs take somewhat lower values than 
those based on the MSTW08 PDFs, particularly 
for the case of Tevatron. This is explained by 
a cumulative effect of the difference in the gluon 
distribution and a s for these two sets. This dis- 
crepancy reduces the statistical significance of the 
constraint on the Higgs mass obtained at the 
Tevatron collider, cf. Refs. [32133) . For the con- 
clusive interpretation of the upcoming LHC re- 
sults it would be useful to reduce the spread of 
the predictions. 

In summary, the ensemble of nucleon PDF sets 
is maintained at the moment and is ready for phe- 



nomenological studies at the LHC. The PDFs are 
tuned to a variety of the hard-scattering data, 
including the most recent ones, and provide an 
accuracy of the LHC standard candle processes 
at the level of several percent [34] . Nonetheless 
the predictions based on the different PDFs are 
not in perfect agreement, presumably due to the 
distinct details in the theoretical treatment of the 
data employed by different groups. Clarification 
of these aspects and consolidation of the predic- 
tions will essentially improve the precision of the 
future LHC findings. 
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